UNCLASSIFIED 


RepA&jLiced. 
by  the 


ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


I 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  puipose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government;  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


u  i  »  < v>  i  it 


63~/~3 


i  286  984 

O 

CD 

00 

AIRPORT  RUNWAY  CAPACITY 
AND  DELAYS  — I 

by 


*  Robert  M.  Oliver 

CD 


OPERATIONS  RESEARCH  CENTER 

INSTITUTE  OF  ENGINEERING  RESEARCH 


RESEARCH  REPORT  25 
28  AUGUST  1962 
I.E.R.  172-28 


i 


UNIVERSITY  OF  C  A  L  I  F  0  R  N  I  A  -  B  E  R  K  E  L  E  Y 


AIRPORT  RUNWAY  CAPACITY  AND  DELAYS  —  I 


by 


Robert  M.  Oliver 
Operations  Research  Center 
University  of  California,  Berkeley 


28  August  1962 


Research  Report  25 


This  research  has  been  partially  supported  by  the  Office  of  Naval  Research 
under  Contract  Nonr-222(83)  with  the  University  of  California.  Reproduc¬ 
tion  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United  States 
Government. 


ABSTRACT 


When  both  landing  and  departing  aircraft  use  a  single  runway 
the  flow  capacity  and  the  delays  to  aircraft  depend  on  the 
operating  rules  which  govern  use  of  the  runway.  In  this  paper 
the  author  analyses  those  situations  where  landing  aircraft 
have  priority  over  departures  but  where  departures  must 
satisfy  several  minimum  spacing  requirements  before  they 
can  be  interposed  between  landings.  By  assuming  that  inter - 
arrival  times  of  landing  aircraft  are  independently  sampled 
random  variables,  from  a  translated  exponential  distribution, 
solutions  are  obtained  for  the  total  movement  capacity  and  the 
delays  to  departures  at  a  single  runway. 


AIRPORT  RUNWAY  CAPACITY  AND  DELAYS _ I 


!•  Introduction 

In  recent  years  demands  for  air  service  have  increased  at  a  rate  which 
often  exceeds  the  rate  at  which  new  facilities  and  new  control  procedures  can 
be  designed  and  installed  to  handle  these  traffic  flows.  In  order  to  provide 
economical  designs  and  efficient  control  procedures,  one  needs  to  examine 
the  delay  and  congestion  problems  which  arise.  This  is  especially  true  when 
terminal  areas  operate  near  saturation,  i.  e.  ,  as  peak  flow  rates  of  landing, 
taxiing  and  departing  aircraft  approach  the  capacity  of  runway  facilities. 

The  theory  of  queues  and  those  of  stochastic  processes  in  general  have 
suggested  many  fruitful  lines  of  approach  to  these  delay  questions.  Since  the 
probability  distributions  and  average  delays  are  extremely  sensitive  to 
changes  in  arrival  rates  when  the  servicing  facility  is  close  to  its  capacity, 
the  need  for  accurate  calculations  of  capacity  is  evident.  Moreover,  this 
theoretical  capacity  reflects  an  upper  bound  on  the  actual  number  of  landing 
and  departing  aircraft  which  can  be  serviced  per  unit  time  over  long  periods 
of  time.  Once  movement  capacities  have  been  obtained  it  is  sometimes 
possible  to  find  delays  to  landing  and  departing  aircraft.  This  paper  discus- 
ses  both  capacity  and  delay  questions. 

Assume  that  landing  aircraft  arrive  at  random  over  a  runway  threshold. 
Under  these  conditions  and  in  the  absence  of  any  departing  aircraft,  the  run¬ 
way  capacity  is  determined  by  the  length  of  time  an  airplane  is  on  the  runway. 

If  the  service  time  of  the  runway  (i.  e.  ,  from  -over  threshold"  to  "off  runway") 
is  a  random  variable  with  mean  value,  A,  the  steady- state  capacity  service 
rate  of  the  runway  is  A*1.  These  statements  are  also  true  when  only  departing 


aircraft  are  allowed  to  use  the  runway.  In  each  case,  the  capacity  or  steady- 
state  flow  rate  through  the  runway  is  the  reciprocal  of  the  mean  service  time. 

Even  if  the  dominant  service  times  occur  on  the  runway  these  simple 
arguments  do  not  extend  to  mixed  operations  for  several  reasons:  the  flow 
stream  of  departing  aircraft  may  not  operate  independently  of  the  stream  of 
landing  aircraft.  Interdeparture  times  may  no  longer  be  independently 
sampled  random  variables  and  the  runway  occupancy  times  of  landing  aircraft 
may  differ  widely  from  those  of  departing  aircraft. 

In  this  paper  we  assume  that  demand  for  service  is  random  in  the  sense 
that  the  interarrival  times  of  landing  aircraft  are  independently  sampled  from 
known  probability  distributions.  This  is  not  so  restrictive  as  the  statement 
that  demand  for  service  is  Poisson;  on  the  other  hand,  it  does  not  allow  one 
to  incorporate  interarrival  times  which  depend  on  the  history  of  earlier 
arrivals. 

For  the  purposes  of  our  analysis,  the  time  scale  t  >  0  can  be  sub¬ 
divided  into  two  types  of  intervals:  blocked  intervals  when  a  departing  air¬ 
craft  is  forbidden  use  of  the  runway,  and  unblocked  intervals  when  he  may 
use  the  runway.  (See  Figure  1.) 

In  Section  2  we  obtain  the  distribution  of  lengths  of  blocked  and  un¬ 
blocked  periods  when  every  point  in  an  unblocked  period  is  removed  by  more 
than  a  constant  separation  X  from  a  preceding  airplane  and  more  than  a 
constant  separation  Y  from  a  following  airplane.  In  Sections  3  and  4 
expressions  for  the  capacity  servicing  rate  of  a  runway  with  both  landings 
and  take-offs  are  developed  from  an  asymptotic  expansion  of  the  average 
count  of  departures  interposed  in  an  interval  t.  Delays  to  departing  air- 
craft  are  studied  in  the  final  sections. 
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LANDINGS  AND  DEPARTURES. 


The  work  of  earlier  authors  is  related  to  the  material  presented  in  this 
paper.  BlumsteuJ  ^  has  discussed  the  operational  capacity  of  a  runway  used 
by  landing  and  departing  aircraft.  In  his  models  the  long-run  or  steady- state 
interposition  rate  of  departing  aircraft  is  found  when  successive  departures 
satisfy  a  Markov  property.  The  renewal  arguments  of  this  paper  allow  more 
general  interposition  rules,  explicit  expressions  for  the  count  of  interpositions 
in  a  finite  interval  of  time  as  well  as  formulas  for  the  delay  distributions  of 
departing  aircraft. 

The  result  in  Section  2  is  a  slight  modification  of  the  ordinary  gap  and 
block  process  in  a  traffic  stream.  The  reader  is  referred  to  the  papers  by 
Tanner^,  Weiss  and  MaraduliJ10^ ,  and  Oliver^5  ^ ,  and  references  contained 
therein  for  more  detailed  treatment  of  this  well-known  problem. 

The  delays  in  the  queue  of  departing  aircraft  is  similar  to  a  paper  by 

[8] 

Tanner  dealing  with  the  interactions  of  two  opposing  Poisson  streams  of 
Interference  and  queues  arise  when  there  is  simultaneous  demand 
for  service  at  a  facility  that  can  serve  only  one  stream  at  a  time.  In  Tanner's 
paper,  priority  is  determined  by  the  first  occupant  of  the  facility.  Recent 
work  by  Gaver^2^  obtains  some  important  results  for  priority  queueing  models 
where  customers  queue  for  gaps  in  a  primary  traffic  stream.  Analogies  be¬ 
tween  his  and  our  model  are  discussed  in  the  text. 

The  notation  used  in  this  paper  is  similar  to  that  used  in  an  earlier 
paper  on  the  distribution  of  blocks  and  unblocked  periods  in  a  traffic  stream?"*^ 
Interevent  density  distributions  are  denoted  by  small  letters.  Capital  letters 
denote  the  probability  that  the  random  variable  is  greater  than  or  equal  to  the 
argument.  Average  values  and  variances  of  the  density  distributions  are  shown 
by  v  and  a2  with  subscripts  referring  to  the  appropriate  density  distribution. 
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X  and  p,  with  or  without  subscripts,  refer  to  steady-state  arrival  and 
service  rates.  Laplace  transforms  are  denoted  by  a  tilde  (  ~  )  over  the 
function  being  transformed  and  6(t)  is  the  dirac  delta  function.  Subscripts 
d  and  i  refer  to  departing  and  landing  aircraft  respectively. 


* 
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2.  Blocked  and  Unblocked  Periods 


In  order  to  study  the  distribution  of  blocked  periods  in  which  departing 

aircraft  are  forbidden  use  of  the  runway,  we  assume  that  aircraft  pass  over 

the  threshold  of  a  runway  at  instants  of  time  rn,  r,,  t,,  r  ...  .  and  that 

u  i  i.  n 

the  intervals  tft  =  Tn+i  -  between  arrivals  are  independently  sampled 
from  a  common  density  distribution  a(t)  with  mean  value  v  and  variance  <r2. 

If  Tq  =  0  and  n  aircraft  pass  over  the  threshold  of  the  runway,  the 
last  one  arriving  at  the  instant  t  ,  an  "arbitrary"  point  has  the  probability 
dt  •  of  landing  in  an  interval  dt.  The  arbitrary  point  must  necessarily 

land  in  an  interval  which  is  bounded  on  the  left-hand  side  by  some  r.  and  on 
the  right  by  some  t.^;  0  ^  j  n-1.  Since  the  arbitrary  point  is  more 
likely  to  be  covered  by  long  than  by  short  interarrival  intervals,  the  probability 
density  distribution,  v(t),  of  the  interarrival  time  which  covers  the  arbitrary 
point  is  proportional  to  t  and  a(t)  as  n  becomes  large^^ : 

v(t)  =  v~  *ta(t)  (1) 

This  distribution  has  mean  and  variance 


(2) 


a2  =  „-V3>  - 


where  i/n^  is  the  n^1  moment  (about  zero)  of  a(t). 

The  conditional  probability  density  distribution  u(r/t)  that  the  spacing 
from  an  arbitrary  point  to  the  following  (or  preceding)  aircraft  is  between  r 
and  r  +  dr  given  that  the  point  is  covered  by  an  interval  of  length  t  is  the 


-6- 


rectangular  distribution  over  (0,t] 


u(r/t)dr  = 


0  <  r  <  t 


(3) 


and  the  marginal  density  distribution  of  spacings  from  the  arbitrary  point  to 
either  a  following  or  preceding  airplane  is  the  " starting -at-random"  or 
"next  arrival"  density*  u(t)  =  i/_1A(t)  which  has  the  following  physical 
interpretation:  If  one  arrives  at  the  runway  threshold  at  an  instant  of  time 
which  is  uncorrelated  with  the  arrival  times  {rn  }  the  probability  density 
distribution  of  wait  to  the  first  airplane  which  passes  over  the  runway  threshold 
is  u(t) •  The  conditional  probability  density  distribution  v(t/r)  that  the 
arbitrary  point  is  covered  by  an  interval  of  length  t  given  the  spacing  r  to 
the  first  following  (or  preceding  vehicle)  is  therefore 


v(t/r )  =  U<r('t>V<t>  =  iW 

v~  *A(r) 


r  <  t 


(4) 


Define  unblocked  points  as  the  collection  of  points  which  are  separated 
by  more  than  an  interval  Y  from  the  first  preceding  airplane  and  more  than 
X  from  the  first  following  airplane.  All  other  points  on  the  real  axis  are 
blocked.  The  probability  that  an  arbitrary  point  will  land  in  an  unblocked 
period  is  the  joint  probability  that  the  spacing  to  the  vehicle  on  the  left  is 
greater  than  or  equal  to  X  and  that  the  spacing  to  the  vehicle  on  the  right  is 

greater  than  or  equal  to  Y.  Hence,  from  the  definition  and  Equation  (4)  we 
get 


Reference  7,  p.  12. 
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Pr  -^arbitrary  point  lands  in  an  unblocked  period 
oo  p  oo 

\  v(r  +  t/t)  u(t)drdt  (5) 

X  Y 

=  U(X  +  Y) 

In  a  traffic  context  this  probability  is  called  the  transparency  of  the  landing 
stream,  i.  e.  ,  the  fraction  of  time  the  landing  stream  can  accept  departures. 
The  probability,  V(t),  that  an  arbitrary  point  is  covered  by  an  interarrival 
time  greater  than  or  equal  to  t,  is  greater  than  U(t),  in  the  former  case, 
we  allow  one  vehicle  to  be  arbitrarily  close  to  the  arbitrary  point  while  we 
do  not  in  the  latter  case.  In  general  we  find  that 

V(t)  -  U(t)  =  tu(t)  >  0 

and  it  is  only  in  the  Poisson  case  (exponentially  distributed  interarrival  times) 
that  the  probability  of  landing  in  an  unblocked  period  is  equal  to  the  product  of 
the  independent  probabilities  that  the  spacing  to  the  preceding  airplane  is 
greater  than  Y  and  to  the  following  is  greater  than  X. 

When  departing  aircraft  must  meet  minimum  spacing  requirements  to 
preceding  and  following  landings  as  well  as  departures,  the  length  of  a  blocked 
period  consists  of  the  sum  of  three  parts:  (i)  an  interval  X  which  precedes 
the  arrival  of  the  first  of  several  landing  aircraft,  each  with  interarrival  time 
less  than  X  +  Y;  (ii)  the  length  of  an  interval  from  this  first  arrival  to  the  first 
arrival  which  begins  a  gap  greater  than  X  +  Y;  and  (iii)  an  interval  Y  which 
follows  the  last  mentioned  arrival.  Each  unblocked  period  must  be  covered  by 
an  interarrival  headway  between  landing  aircraft  which  is  greater  than  or  equal 
to  Z  =  X  +  Y.  Clearly,  the  distribution  of  lengths  of  unblocked  periods  does  not 
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depend  on  the  location  of  the  unblocked  period  between  two  arrivals  but  only 
on  the  distribution  of  intervals  which  exceed  Z.  The  density  distribution  of 
lengths  of  unblocked  periods,  h(t),  has  a  transforrJ5^ 

h(s)  =  eS  A'  <Z)[a(s)  -  a(s;Z>]  (6) 

where  a(s)  is  the  transform  of  a(t)  and  S(s;Z)  =  f  *  a<t)e- Stdt.  We  denote 

J0 

by  f(t)  the  density  distribution  of  lengths  of  blocked  periods,  i.  e.  ,  those 
points  which  lie  outside  the  minimum  spacing  requirements  to  adjacent  air¬ 
craft.  It  was  also  shown  that  this  density  distribution  has  a  Laplace  transform 


~f<s)  =  -iLlZA(z) 
1  -  a(s;Z) 


(7) 


For  future  reference  we  also  note  that  the  first  and  second  moment  of  f(t)  can 
be  written  in  the  form 


vi  =  vK  ^ZMl  -  U(Z)  ) 


(8) 


:^-2t/fZ  +  Z2  +  A-1(Z)|I/(2)  -  J  t2a(t)dt  j 

J 


(9) 
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3.  Interpositions 


It  should  be  obvious  from  a  practical  as  well  as  a  theoretical  viewpoint 
that  the  number  of  departing  aircraft  which  can  merge  into  the  unblocked 
periods  of  the  stream  of  landing  aircraft,  is  not  equal  to  the  difference  be¬ 
tween  a  capacity  service  rate  of  the  runway  and  the  observed  flow  rate  of 
landing  aircraft.  While  this  computation  may  serve  as  an  approximation  in 
certain  situations,  the  simple  analysis  breaks  down  for  reasons  we  have 
already  mentioned.  Departing  aircraft  may  only  be  able  to  merge  with  the 
stream  of  landing  aircraft  when  certain  minimum  separations  are  maintained 
between  adjacent  landing  as  well  as  departing  aircraft. 

To  obtain  the  counting  distribution  of  departing  aircraft  which  can  be 
interposed  between  landing  aircraft,  we  consider  the  hypothetical  situation 
where  an  infinite  queue  of  departing  aircraft  waits  at  the  side  of  the  runway. 
When  an  unblocked  period  appears,  a  departure  is  interposed,  ie.,  is  allowed 
to  take  off.  If  the  unblocked  period  is  long  enough,  there  may  be  several  inter¬ 
positions.  In  no  case,  however,  are  departing  aircraft  allowed  to  come 
arbitrarily  close  to  one  another. 

Interpositions  can  be  counted  in  several  ways.  For  example,  one  can 
specify  intervals  which  begin,  (i)  just  after  a  landing  has  crossed  the  runway 
threshold,  (ii)  just  after  a  departure  is  released,  (iii)  at  random,  (iv)  at  the 
beginning  of  an  unblocked  period,  or  (v)  at  the  beginning  of  a  blocked  period. 

The  count  depends  not  only  on  the  location  of  this  interval,  but  also  on  the 
rules  which  govern  the  precise  instants  when  departing  aircraft  are  released 
into  the  primary  stream. 

In  the  remainder  of  this  paper  we  focus  attention  on  an  interval  (0,  t) 
which  begins  with  either  blocked  or  unblocked  periods.  If  we  start  to  count 
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at  the  instant  when  a  blocked  period  begins,  the  first  departure  is  released 
at  the  instant  that  the  blocked  period  ends,  i.e.  ,  the  unblocked  period  begins. 
No  departures  will  be  released  in  an  interval  (0,  t)  beginning  with  a  blocked 
period  if  the  blocked  period  is  greater  than  t.  If  we  let 


pn(t)  =  Pr 


Jn  interpositions  in  [0,  t)  starting  at 
^the  beginning  of  a  blocked  period 


(10) 


we  have  the  special  case, 


f ( x)  dx  =  F(t) 


(11) 


If  the  first  blocked  period  ends  at  an  instant  r  such  that  r  <  t,  the  count 
in  (0,  t)  equals  the  count  in  [r,t).  But  the  count  in  [r,t)  is  the  count  one 
would  observe  in  an  interval  [0,  t-r)  whose  origin  was  located  at  the  beginning 
of  an  unblocked  period-  Defining  qn(t)  as  the  discrete  counting  distribution 
of  departures  when  the  counting  origin  starts  at  the  beginning  of  an  unblocked 
period,  we  have 


Pn(0  = 


n  >  1 


(12) 


In  order  to  obtain  a  second  renewal  equation  expressing  q^ft)  in  terms  of 

p  (t)  we  consider  the  events  which  lead  to  a  count  of  n.  if  the  origin  is  located 
n 

at  the  beginning  of  an  unblocked  period.  To  do  this  we  use  still  another 
conditional  probability  distribution.  If  an  unblocked  period  is  of  length  r,  we 
define 


,  .  I  .  „  fm  interpositions  lunblockedl 

h(m|  r)  =  Pr  jinterval  o£  length  r  j 


(13) 
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and  assume  that  h(mj  r)  is  1  if  operating  rules  allow  m  interpositions  in 
r  or  0  if  they  do  not.  A  total  count  of  n  interpositions  in  the  interval 
[0,t)  starting  with  an  unblocked  period  of  length  r  can  be  observed  if  m  <  n 
are  observed  in  [0,  r)  and  (n-m)  are  observed  in  [r,t).  A  count  of  n  will 
also  be  observed  in  t  if  r  >  t  and  n  departing  aircraft  can  be  interposed 
in  [0,t).  The  probability  of  n  interpositions,  is  the  sum,  overall  m,  of 
the  joint  probabilities  discussed  above.  Hence, 

t  n 

^  h(m|r)h(r)pn_m(t-r)dr  +h(n|t)H(t)  .  n>  1  (14) 

m=l 


The  solutions  of  qn(t)  (or  Pn(t)  )  can  be  written  as  a  generating  function  in 
terms  of 


h(z;s)  =  ^  J  h(n|t)h(t)e"8tzndt 
n=l  0 


(15a) 


H(z;s)=^T  J  h(n|t)H(t)e'Stzndt  . 
n=l  0 


(15b) 


Multiplying  both  sides  of  Equation  (14)  by  zn  and  summing  over  n  gives 


Q(z;s)  =  y  J  qn(t)e'stzndt  = 
n=l  o 


-1 


h(z;  s)(l-f  ( s)  )  +  sH(z;  s) 


s -  sf (s)h(z;s) 


(16) 


The  analagous  transform  for  P(z;s)  is  obtained  by  multiplying  Q(z;s)  by 
f  (s)  and  adding  Pq(s). 
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In  the  remainder  of  this  paper  we  consider  capacity  interposition  rates 


and  queueing  effects  under  the  specific  assumption  that  a  headway  between 
two  landings  exceeds  t  has  probability 


A(t)  =  1 


0  <  t  <  A 


-  e 


Ai 


(17) 


When  is  the  average  arrival  rate  of  landing  aircraft  and  A^  is  the 
minimum  headway  between  landings,  this  distribution  has  been  found  to  be 

l4l 

a  reasonably  accurate  description  of  aircraft  movements  over  threshold.  1 
We  observe,  in  passing,  that  A^  =  0  corresponds  to  the  case  of  Poisson 
arrivals. 

In  order  to  construct  the  blocked  and  unblocked  periods  in  the  stream 
of  landing  aircraft  we  now  specify  three  constants  A^,  A^  and  A^ :  A 

departing  airplane  on  the  runway  must  at  all  times  maintain  a  headway 
gre  ater  than  or  equal  to  A^  ahead  of  a  following  landing  and  a  headway 
greater  than  or  equal  to  A^  behind  an  earlier  landing.  In  addition,  two 
departures  located  within  a  common  unblocked  period  must  maintain  at 
least  the  headway  A^  between  each  other.  If  we  now  assume  (Figure  1) 
that  the  minimum  size  blocked  period  must  equal  the  sum  of  these  three 
constants,  algebraic  simplifications  in  our  formulas  are  obtained  by  setting 


A  =  AJ4  +  A, .  +  A 


di 


id 


It  may  happen  that  A.  >  A  in  which  case  there  can  be  one  or  more  inter- 


positions  between  every  landing;  we  assume,  however,  that  A  >  A^  .  Since 

A  >  Ad  il  is  also  obvious  that  two  departures  located  in  adjacent  unblocked 
periods  do  not  interfere  with  one  another's  movements. 

Substituting  the  transforms 


'  8  Ai  /  V  1 

0  =  V  (  +  s(l-X.iAi)J 


a(s;A)  +  s(l-\/Ai)^e  S  l- 


sA.  -  sA  - 


^(A-A^ 


into  Equation  (6)  leads  to  the  exponential  distribution. 


h(t)  = 


exp  -  (\tt/  1-X^) 


for  the  length  of  unblocked  periods.  If  the  first  departure  is  interposed  at 
the  beginning  of  an  unblocked  period  and  an  interposition  is  made  every  A^ 
time  units  thereafter  within  the  same  unblocked  period,  the  transform 
h(z;s)  in  Equation  (15a)  can  be  written 


x  ,  ,  ~  zX.  T  -(V+s)A  j-  -(X.+s)A,-| 

\H(z;s)  =h(z;s)  l-e  d  1  -  ze  d 


where  \  =  1  -  X.^)'  >  \  . 


The  Laplace  transform  of  the  average  count  of  interpositions, 


_  pOO 

M(s)  =  \  M(t)e~Stdt 


j*  I  nqn(t)e-Stdt  , 


0  n=0 
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z 


equals  the  partial  derivative  of  Q(z;s)  in  Equation  (16),  with  respect  to 
at  z  =  1.  Substituting  Equation  (6),  (18)  and  (20)  into  Equation  (16)  leads, 
asymptotically,  to  the  expression  for  the  capacity  interposition  rate  of 
departing  aircraft: 


=  Lim 
t  “•‘00 


r  M(t)  = 


Lim 
s  — •*() 


s2M(s) 


(22) 


=  V 


1  r  "i- 1 

exp  -  (X^A-  Aj)  /  1  -  X^)  1  -  exp  -  (X^  /  1  -  X^) 


Since  the  partial  derivative  of  h(z;s)  with  respect  to  z  at  z  =  1,  s  =  0 
equals  the  average  number  of  interpositions  per  unblocked  period  and  the 
average  flow  rate  of  unblocked  periods  equals  pA(A),  jj.  can  be  written  as 
the  product  of  such  terms,  an  observation  which  seems  intuitive  in  retrospect 
but  difficult  to  prove  in  individual  cases.  However,  these  results  hold 
asymptotically  for  a  broad  class  of  statistical  interposition  rules  described 
by  Jewell. 

For  small  X^,  the  capacity  interposition  rate, 

nd  =  [UXi<*-V2)J  +  .  (23) 

is  equal  to  the  capacity  service  rate  for  a  "departures -only"  runway  operation 
minus  a  small  term  which  is  linear  in  X^ ,  the  average  arrival  rate  of  landings 
Figure  2  is  a  plot  of  pd  for  several  values  of  A^,  Ad  and  A  which  are 
encountered  in  typical  runway  operations  with  jet  and  propellor  powered  air¬ 
craft  in  the  U.  S.A. 

A  similar  transform  can  be  obtained  for  the  counting  distribution  p  (t). 
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Landing  Rate,  X, 

Fig.2 -CAPACITY  INTERPOSITION  RATE. 


5.  Delays  to  Departing  Aircraft 


Once  the  arrival  pattern  of  landing  aircraft  and  the  capacity  inter¬ 
position  rate  of  departing  aircraft  is  known,  it  is  possible  to  obtain  an 
expression  for  the  distribution  of  delays  to  departing  aircraft.  In  particular, 
we  are  interested  in  the  distribution  of  elapsed  time  from  arrival  to  inter¬ 
position.  With  some  differences,  which  we  shall  discuss  shortly,  these 
types  of  delay  problems  can  be  studied  under  the  heading  of  priority  queueing 
problems  where  the  lower  priority  (departing  aircraft)  must  wait  for  unblocked 
periods  to  appear  in  the  high  priority  stream  before  service  can  begin. 

If  we  focus  attention  on  those  instants  of  time  when  a  departure  begins 

use  of  the  runway,  the  headway  to  the  earliest  succeeding  departure  (in  the 

absence  of  landings)  is  A^.  Since  the  length  of  unblocked  periods  (Equation 

(19)  )  is  exponential  with  constant  X.  =  V  X-XA  >  ,  the  probability  that 

the  headway  from  the  beginning  of  service  to  the  beginning  of  the  blocked 

-\Ad 

period  exceeds  A^  is  e  Hence  this  is  the  probability  that  another 

departure  precedes  a  landing,  i.e.,  the  service  time  of  the  second  departure 

is  also  A,. 

a 

On  the  other  hand,  the  probability  that  at  least  one  landing  precedes 

‘XAd 

another  departure  is  1-e  .  In  this  case  the  headway  to  the  first  following 

departure  is 

S  =  (SQ)  +  (Sj  +  S2  +...)+  (A)  (24) 

where  Sq  is  the  headway  (  <A^)  to  the  beginning  of  the  unblocked  period, 

S^(i  >  1)  are  those  headways  between  landing  aircraft  which  lie  within  a 
blocked  period  and  the  last  term  is  the  sum  of  the  minimum  clearances  be¬ 
tween  adjacent  aircraft. 


The  probability  density  of  Sq  is 

Pr  {S0  -  *  }  =  ~^T 

1-e 


0  <  t  <  A , 

—  Q 


(25) 


and  since  the  arrangement  of  remaining  terms  in  Equation  (24)  is  not  important, 
the  probability  distribution  of  S  -  SQ  is  the  distribution  of  lengths  of  blocked 
periods.  The  density  distribution  d(t)  of  the  headway  between  successive 
departures  is  now  obtained  by  removing  the  condition  on  the  length  of  SQ  and 
adding  the  probability  that  the  service  time  is  exactly  A^: 

~X  A .  p  d  . 

d(t)  =  e  6(t-Ad)  +  \  Xe"*rf(t-r)dt  (26) 

Jo 

We  again  point  out  that  the  blocked  periods  are  constructed  from  a  landing 
stream  with  interarrival  headways  given  by  Equation  (17).  Substituting 
Equation  (18)  into  Equation  (7)  for  f(s)  and  collecting  terms,  we  obtain  the 
Laplace  transform  of  the  probability  that  the  service  time  of  a  departure 
exceeds  t  as. 


00 


D(s)=J  D(t)e"stdt  =  ^  +  s-\e 


-sA, 


-(\  +  s)Ad 


1 


X  +  s  -  Xe 


+  X« 


-(X+s)A  XA. 


(27) 


As  s  r-*-  0  we  obtain  the  mean  service  time  of  departures 


v 


d 


Lim 
s  — *•  0 


D(s)  =  X’ 1 


X(A-A). 

0 


e'XAd) 


(28) 


which  is  just  the  reciprocal  of  the  capacity  service  rate  in  Equation  (22). 
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The  expected  delays  to  departing  aircraft  can  be  obtained  by  direct 
renewal  arguments  when  we  assume  that  arrivals  to  the  departure  queue  are 
Poisson  with  mean  flow  rate  and  that  first-come,  first-serve  rules  are 

observed  for  departures.  Let  Wq  be  the  probability  that  there  are  no 
departures  on  the  runway  or  waiting  to  be  interposed  between  departures. 
Since,  for  long  run  equilibrium  conditions,  the  expected  number  of  departures 
demanding  service  must  equal  the  expected  number  interposed  in  the  landing 
stream,  we  have 

^dl  =  1  _7ro>t  +  ^ 

Hence  as  t  — ►  oo  (29) 

*o  =  1  -  Xdvd  • 

A  departure  arriving  at  an  arbitrary  instant  of  time  waits  for  use  of  the  run¬ 
way;  this  random  wait  has  a  density  distribution  w(t).  If  the  arrival  finds 
one  or  more  departures  waiting,  the  probability  density  of  his  wait  equaling 
t  is 

©  ^ 

^  D(r)w(t-r)dr  t  >  0  (30) 

since  p^D(r)  is  the  probability  density  that  the  departure  at  the  head  of  the 
queue  must  wait  a  time  r  for  use  of  the  runway.  If  a  departure  arrives 
when  no  other  departures  are  waiting,  he  may  occupy  the  runway  at  once, 
i.  e.  ,  the  wait  is  zero,  or  he  may  have  to  wait  until  the  current  blocked 
period,  of  minimum  size  A,  passes.  In  the  former  case,  the  probability  of 
landing  in  an  unblocked  period  is,  by  a  long  run  argument,  equal  to  the 
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fraction  of  time  occupied  by  unblocked  periods,  i.  e. ,  the  transparency  of 
the  stream. 


U(A)  =  (1+X.i^)'1 


(31) 


while  the  complementary  probability  is  the  chance  of  landing  in  a  blocked 
period.  If  the  departure  arrives  during  a  blocked  period  the  probability 
density  of  headways  to  the  end  of  the  blocked  period  is  v~  1F(t) .  Adding  up 
the  probabilities  of  exclusive  events  which  lead  to  a  wait  of  t,  we  finally 
get  the  density  distribution  of  wait  for  use  of  the  runway: 


w(t) 


uxdvd  r 

-  TT\u{  [' 


6(t)  +  XF(t) 


D(r)w(t-r)dr  . 


(32) 


The  Laplace  transform  of  Equation  (32), 


.  l+XF(s) 
w(s)  =  - 

1  -XdD(s) 

is  a  modification  of  the  Pollaczek-Khintchine  formula^ which  accounts  for 
the  fraction  of  time  when  departures  arrive  too  close  to  landings  for  safe 
take-off  operations. 

[2] 

Gaver  ,  has  obtained  a  similar  result  for  the  waiting  time 

distribution  of  vehicles  which  queue  for  gaps  between  blocked  periods.  His 

blocked  periods  would  have  (in  our  runway  model)  a  minimum  size 

A-  Ad  =  Adi  +  Aid  ;  lying  between  each  blocked  period  is  a  gap  which  can  be 

used  by  one  (or  more)  departing  aircraft  provided  its  size  exceeds  A,,  the 

a 

minimum  headway  required  by  departures.  As  a  pre-emptive  priority  repeat 


uVd 

1  +  \vr 


(33) 
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queueing  model,  he  finds  that  a  departure's  "completion-time"  consists  of 
an  alternating  sequence  of  pre-emptions  by  landings  (blocked  periods)  and 
insufficient  gaps  (unblocked  periods  <  Ad)  before  it  is  actually  allowed  onto 
the  runway. 

Let  b(s)  be  the  transform  of  the  blocked  period  distribution  obtained 
when  interlanding  times  are  given  by  Equation  (17)  and  we  substitute 
A^  +  Ajd  for  Z  in  Equation  (6)  and  (7).  By  using  the  identity 


f  (s)  = 


-(s+\)A 
e  b(s) 


1 


\b( 


>(b) 

+  s 


-(UsT^ 


it  is  possible  to  obtain  Gaver's  result, 


- s  A,  _  -1 

w(s)e  [  d(s)] 


1  -  Vd 

l-XdD(s) 


l+XB(s) 

1+X,/b 


(34) 


(35) 


for  the  transform  of  the  waiting  time  that  elapses  from  the  instant  a  departure 
joins  the  queue  until  it  is  at  the  head  of  the  line. 

The  expected  delay,  wd>  in  the  departure  queue  can  be  obtained  by 
direct  integration  of  Equation  (32): 


wd  * 


X,„J 2) 
d  d 


\v 


(2) 


2t1-Xdvd)  2(1+Xvf) 


(36) 


As  we  have  already  mentioned,  the  general  characteristics  of  the  Pollaczek- 
Khintchine  formula  are  apparent  in  the  first  term.  The  second  term  is  the 

*  Reference  [2],  Section  7. 
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additional  contribution  of  blocked  periods  to  average  delay  of  departures. 
One  might  be  tempted  to  conclude  that  the  effect  of  the  blocked  periods  is 
small, since  the  major  contribution  to  average  delay  (with  high  departure 
flow  rates)  is  due  to  the  left  hand  term  in  (36)  involving  only  the  first  and 
second  moment  of  the  service  time  distribution.  However,  it  must  be 
remembered  that  both  v ^  and  i/j  ^  depend  upon  the  length  of  blocked 
periods.  (Equations  (26)  and  (27)). 

Numerical  calculations  can  now  be  completed  by  substituting  the 
moments  of  the  service  time  distribution  in  Equation  (26)  and  moments  of 
blocked  period  lengths  which  can  be  obtained  by  substituting  Equation  (17) 
into  Equations  (8)  and  (9)-  Since  the  algebraic  expressions  are  involved, 
we  now  turn  to  low  flow  rate  approximations  and  the  special  case  of  Poisson 
landings  and  departures. 


-22- 


6.  Low  Flow  Approximations 


As  either  the  arrival  rate  of  landing  or  departing  aircraft  becomes 
small,  the  expressions  for  capacity  servicing  rates  and  delays  simplify 
considerably.  For  low  arrival  rates  of  landing  aircraft,  i/^  — ►  and 

V d^  —  A^  .  The  transform  of  the  waiting  distribution  (Equation  (33) ) 
approaches  that  of  a  constant  service  queue  since  the  restriction  on 
immediate  interposition  of  an  arbitrary  departure  is  due  to  the  minimum 
separation  A^  between  departures  but  not  to  the  higher  priority  landing 
aircraft.  In  particular,  the  average  delay  to  departing  aircraft  becomes 


Liim  w 


XdAd 


2(l-XdAd) 


(37) 


On  the  other  hand,  as  the  average  arrival  rate  of  departures  fall  off 
to  zero,  departure  delays  are  primarily  determined  by  the  presence  or 
absence  of  landing  aircraft  while  the  effect  of  interference  between  depar¬ 
tures  is  small.  Expansion  of  Equation  (33)  gives 


w(s)  = 


(38) 


-MA-  Aj) 


-sA. 

(X+s)lX+s-Xe  1 

na; -is  ^.ta=a;j 


) 


s(l+XA^)(X  +  s -Xe  +  Xe  e  J 


d  +  0(Xd) 


The  term  in  square  brackets  is  the  transform  of  the  wait  for  an  unblocked 

period.  The  limiting  value  of  the  average  delay  to  departing  aircraft  is  there- 

$ 

fore  equal  to  the  average  wait  for  an  unblocked  period  in  the  landing  stream. 

jjC  __ 

Equation  (42),  p.207,  and  Equation  (47b),  p.208.  Reference  [5]  . 
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When  landings  are  Poisson,  A^  =  0  and  X.  =  X^  in  Equation  (17).  If 
we  also  assume  that  and  A^  are  small  compared  to  A^,  the  blocked 

periods  and  the  time  between  interposition  of  departures  have  identical 
distributions  with  average  value  and  second  moment  equal  to 


Substituting  these  expressions  into  Equation  (36)  gives 


\i  +  \d 
Wd  =  '  ^ 


XiAd  . 

e  -  y^d -1 

,  ,  v  V  -1  V^d  J 

1  4*  XA,  -  e 


(40) 


for  the  average  delay  to  departures.  For  positive  values  of  X^  less  than 
the  capacity  service  rate,  the  term  in  square  brackets  is  also  positive  and 
a  plot  of  w^  as  a  function  of  X^  has  the  characteristic  shape  of  the  delay 
to  arrivals  in  a  single  channel  queue. 
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7.  Summary 


The  analysis  of  aircraft  service  rates  and  delays  has  focused  on  mixed 
runway  operations  where  low  priority  departures  are  interposed  between 
high  priority  landing  aircraft  and  where  minimum  separation  standards  are 
required  between  all  aircraft.  Capacity  interposition  rates  of  the  servicing 
system  have  been  found,  since  they  yield  an  upper  bound  to  the  interposition 
or  total  movement  rate  which  can  be  sustained  over  long  periods  of  time. 
Even  in  the  absence  of  priority  queueing  models  which  incorporate  specific 
holding  patterns  for  landings  or  departures,  one  can  expect  to  find  measures 
of  cost  and  delay  built  around  such  capacity  figures.  In  addition,  they  can 
also  be  used  to  assess  the  trade-off  between  landing  and  departure  service 
rates. 

Actual  operations  on  a  runway  are  complicated  by  several  factors 
which  we  did  not  include.  In  the  first  place,  the  concave  (convex)  time- 
distance  trajectories  of  landings  (departures)  may  begin  or  terminate  at 
several  rather  than  one  point  due  to  the  introduction  of  multiple  taxiway s. 

The  shape  of  these  curves  will  also  depend  on  the  type  of  taxiways,  landing 
speeds  and  accellerating  and  braking  characteristics  of  aircraft. 

In  general,  the  time  (or  distance)  an  airplane  occupies  the  runway  will 
be  a  random  variable.  As  the  real-time  control  of  runway  operations  be¬ 
comes  more  precise  through  the  use  of  semi-automatic  surveillance  and 
computing  equipment,  minimum  separation  standards  may  depend  upon  air¬ 
craft  type  within  a  priority  class  rather  than  be  constant  as  we  have  assumed. 
In  addition,  simultaneous  minimum  distance  and  headway  criteria  may  be 
imposed. 

We  have  also  assumed  that  delays  to  congestion  are  caused  by  runway 
operations.  As  aircraft  flight  characteristics  and  instrument  flight  rules 
change,  dominant  service  operations  may  appear  in  the  glide  path  or  along 
the  airways  which  lead  into  a  terminal  area  rather  than  on  the  ground.  How¬ 
ever,  one  believes  intuitively  that  many  of  the  results  obtained  in  these 
simpler  models  will  extend  to  those  which  include  the  effects  mentioned 
above. 
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